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Abstract 



The virtual corrections to the production cross section of a Standard Model Higgs 
boson are computed up to order ot\. Using an effective Lagrangian for the limit 
Mt — * oo, we evaluate the relevant massless two-loop vertex diagrams by mapping 
them onto three-loop two-point functions, following a method recently introduced by 
Baikov and Smirnov As a result, we find a gauge-invariant contribution to the 
| total Higgs production cross section at NNLO. 

o 

^ ! 1 Introduction 

Q_i| The Standard Model of Elementary Particle Physics (SM) has been verified in many of 

its details with enormous precision over the last 20 years. However, several fundamental 
questions remain unanswered. Perhaps the most important unresolved problem concerns 
the origin of the particle masses. 

In the SM, the underlying mechanism for mass generation is spontaneous breakdown of the 
electro-weak gauge symmetry. In its minimal version (which we will assume throughout 
this paper) it predicts a single, as yet undetected physical particle, the Higgs boson. It 
is determined to be electrically neutral and of spin zero. Its mass, however, is a free 
parameter of the theory. 

The extensive searches for the Higgs boson at particle colliders have set a lower bound 
on its mass, the latest update yielding a limit of around 108 GeV ||. On the other 
hand, theoretical predictions for physical observables depend on the Higgs mass through 
radiative corrections. Comparison of the existing experimental precision data with their 
theoretically predicted values allows to derive a most probable range for the Higgs mass 
which turns out to be roughly between 100 and 200 GeV • 

With LEP being close to its maximum possible energy, the attention concerning Higgs 
search is turning towards future experiments at hadron colliders, in particular LHC, sched- 
uled for the year 2007, or already Tevatron's Run II, starting in 2001. 
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Figure 1: Leading order diagram to the process gg — > H: (a) in the full and (b) in 
the effective theory. The "(8)" denotes the effective vertex of Eq. (pi) . 



The dominant production mechanism for a Higgs boson with a mass below 1 TeV at the 
LHC will be through gluon-gluon fusion (for a review see ||). The coupling of the gluons 
to the Higgs boson is mediated through a quark loop, Fig. ffl(a). In the heavy quark limit, 
the corresponding form factor becomes independent of the quark mass. Thus, this process 
can be used, for example, to count the number of heavy quarks that may exist beyond the 
third generation. 

The current theoretical prediction for this reaction carries an uncertainty of about a factor 
of 1.5 to 2. It is therefore important to improve on the theoretical accuracy. In this paper 
we provide a gauge invariant ingredient to the complete next-to-next-to-leading order 
prediction, namely the virtual corrections up to order a*. The calculation is, to our 
knowledge, the first application of a recently introduced method that allows to relate the 
relevant set of vertex diagrams to the more familiar class of three-loop two-point functions. 



2 Effective Lagrangian 

As it was mentioned before, the coupling of the gluons to the Higgs boson is mediated 
through a quark loop, Fig. ||(a). Since all quarks except for the top are much lighter than 
the current lower limit on the Higgs mass, we will neglect their masses in the following. In 
this case, the top quark is the only one that couples directly to the Higgs boson, because 
the Higgs-fermion vertex is proportional to the fermion mass. The leading order result 
has been known for quite a while M]. At the parton level it reads: 
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where s is the partonic cms energy and G F is the Fermi coupling constant. a s is the strong 
coupling constant which depends on the renormalization scale /i. Mt is the pole mass of 
the top quark, and Mh is the Higgs mass. In order to arrive at the cross section for hadron 
collisions, a hQ has to be folded with the gluon distribution functions. 
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The currently favored values for the Higgs mass appear to be not too different from the 
top quark mass. However, since the threshold for open top production is at Mh = 2M t , 
an expansion in terms of small Higgs mass is expected to work well for Mh < 2M t . In 
fact, for the most interesting mass range of 100GeV< Mh <200GeV it appears that at 
next-to-leading order in a s the complete result is excellently approximated by the leading 
term in an expansion in M^/M 2 [|],||]. Therefore we think it is reasonable to adopt this 
limit also at next-to-next-to-leading order (NNLO). 

All calculations in this paper have been performed using dimensional regularization in 
D = 4 — 2e space time dimensions. Unless stated otherwise, renormalized expressions are 
to be understood in the MS scheme, bare ones will be marked by the superscript "B". 
Furthermore, all the quantities used in the following refer to the five flavor effective the- 
ory. For example, by a s we mean the running coupling constant with five active flavors, 

ai 5 V). 

The most convenient way to obtain the leading term in M^/M 2 is to use the follow- 
ing effective Lagrangian for the Higgs-gluon interaction, where the top quark has been 
integrated out: 
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Here, G^ is the gluonic field strength tensor in the effective five flavor theory, and v 
is the vacuum expectation value of the Higgs field, related to the Fermi constant by 
v = (y^Gp)" 1 / 2 . The coefficient function C\ has been computed in @ up to 0{a A s ). In 
order to obtain the cross section for gg — » H with NNLO accuracy, however, it will only 
be needed up to O(of) @: 
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where It = ln(^ 2 /M 2 ), with Mt the on-shell top quark mass. Here and in the following, 
the number of (light) flavors n\ will eventually be set to five. 



The renormalized operator in Eq. (|2|) is related to the bare one through 

1 
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where 



M = exp [e (-7 B + In Att)} , 
and (3{a s ) governs the running of the strong coupling constant: 

/i 2 ^2«s = a s 0{a s ) . 
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Its perturbative expansion is known up to O(q^) [|T(|, but for the present purpose the 
terms up to C(a 2 ) are sufficient: 
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Using the Lagrangian of Eq. (0) for the Higgs-gluon interaction instead of the full SM, 
the number of loops reduces by one. For example, in leading order one obtains the tree 
diagram shown in Fig. 0(b). The corresponding expression for the partonic cross section 
is 

and coincides with the limit r — > oo of Eq. ([]]), of course. 

Going to higher order, one has to compute virtual corrections to diagram Fig.[l](b). How- 
ever, these will contain infra-red and collinear divergences. The infra-red divergences will 
be canceled when adding the corrections from real radiation of quarks and gluons, but the 
sum will still have collinear divergences. The latter will disappear only when the Altarelli- 
Parisi splitting functions are taken into account up to the appropriate order fTlfl . At the 
next-to-leading order, a full calculation has been carried out in ||. The correction terms 
increase the cross section by about a factor of 1.5 to 2 in the relevant Higgs mass range. 



At NNLO, only a few ingredients for the full answer are available: In [12] the one-loop 



amplitude for the radiation of a single quark or gluon in gluon-gluon, gluon-quark, and 



quark-quark fusion was obtained, and [13] contains the tree- level amplitude for the double- 
emission of gluons and quarks in these reactions. Both of these results still have to be 
integrated over the corresponding two- and three-particle phase space, which certainly is 
a rather non-trivial task by itself. 

In this paper, we want to add the virtual two- loop corrections to the list of available 
knowledge at NNLO. Together with the phase-space integrated expressions for the real 



radiation and the Altarelli-Parisi splitting functions [11], one will then be able to arrive 
at a finite result for the partonic cross section. In order to obtain a physically accessible 
quantity, one needs to fold this partonic result with the corresponding parton distribution 
functions. However, they have yet to be evaluated to the appropriate order in a s (see, 
e.g. EJ" 



3 Calculation of the two-loop diagrams 

Examples of diagrams contributing to the virtual two-loop corrections to the process gg — > 
H are shown in Fig. ^. The Higgs boson couples to the effective vertex resulting from the 
Lagrangian of Eq. (^). The two gluons are on-shell (p\ = p 2 , = 0), which is why — after 
extraction of the tensor structure — the diagrams depend only on one kinematic variable, 
(Pi + P2? = q 2 = 

In addition to their topology and the power of the denominators, the resulting integrals 
can be classified by the power of irreducible numerators, i.e. invariants of momenta that 
cannot be expressed in terms of denominator s.[] In [15, 16|, the integrals with unit (or 



zero) power of denominators and low powers of irreducible numerators were evaluated 
using Feynman parameterization and dispersion techniques. 



In 1 17] recurrence relations based on the integration- by-parts algorithm [18] for the pla- 



nar two-loop integrals were derived, reducing them to convolutions of one-loop integrals. 



lr There is a freedom in choosing the specific invariants, of course. 
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Figure 2: Sample two-loop diagrams contributing to gg — > H at NNLO. The right 
vertex stands for the effective coupling of the Higgs particle to gluons. The bubble 
in diagram (e) represents any quark except for the top quark. 



These relations allow to compute any such planar diagram with arbitrary powers of the 
denominators and irreducible numerators. 

In our case, however, we also need to compute non-planar diagrams, e.g. Fig. ^ (b). For this 
reason, we follow an algorithm that has recently been published by Baikov and Smirnov [|l| . 
It relates the recurrence relations for /-loop integrals with n + \ external legs to the ones 
for (I + l)-loop integrals with n external legs. Here we have n = I = 2, and thus the 
massless two-loop vertex diagrams of Fig. ^ are mapped onto massless three- loop two- 
point functions. The algorithm to compute the latter ones is known [|ll| and implemented 
in the computer program MINCER [l9|| , written in FORM [pQ] . Following the recipe of |M, 
we modified the MINCER routines such that they are applicable to the class of two-loop 
three-point functions at hand. For the generation of the diagrams we used QGRAF |H] as 
integrated in the program package GEFICOM [22|]. 



The only integral that can not be reduced to convolutions of one-loop integrals in this 
approach is the non-planar one with all propagators appearing in single power, and with 
numerator equal to one. However, the result for this integral is known as an expansion in 
e up to its finite part [15]. 



As a check of our setup we re-did the calculation of the electro-magnetic quark form factor 



in QCD to two loops and found full agreement with [16]. We performed this calculation in 
a general gauge and explicitly checked its gauge parameter independence in this way. 
We also computed the two-loop three-gluon vertex in R^ gauge with two gluons on-shell 



and found agreement with |17,23,24|. Finally, the calculation of the present paper was also 
performed in R^ gauge and we verified that the gauge parameter dependence disappears 
in the sum of all diagrams. 



2 I acknowledge the kind permission by the authors of GEFICOM to use this program. 
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4 Results 



The virtual cross section for the process gg — > H can be written as 

*" - 4 ^ * (1 - *> ( l-fl/k.)/. ) 256(1^6? £ • < 9 ) 

where 

X = e ^(pi,A 1 )e^(p 2 ,A 2 )^(Pi,P2), (10) 

and Ci is the coefficient function given in Eq. (Q). The factor 1/256/(1 — e) 2 comes from 
the average over initial state polarizations and color. The kinematical constraint on the 
momenta is 

Pi+P2 = q, Q 2 = M^, p\=pl = Q. (11) 

£ a '^(p, A) is the polarization vector of a gluon with momentum p and polarization A, and 
a, b, respectively fj,, v are the SU(3)-color and the Lorentz indices. Adopting a covariant 
R% gauge, the polarization vectors obey the relation 

Y j e a ;(p,X)e b u ( P ,X) = -9^5 ab . (12) 

A 



The tensor A°^, may be written as 
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Af v = l a • V?)9iit> + bpivp 2f , + cp ltl p 2 u + dp^piv + ep 2 pP2v) • (13) 

H 

The last two terms seem to violate gauge invariance. However, to lowest order they vanish, 
and starting from next-to-leading order, their contribution to the squared matrix element 
gets canceled by the diagrams with ghosts in the initial state. We explicitly checked this 
cancelation by computing these ghost diagrams to the corresponding order. Further, we 
have 

a = -b, (14) 

and thus the c term in (|D|) does not contribute to the total rate. Therefore we will only 
quote the result for a in the following. Its loop expansion can be written as 

a= i + < a (D + f<y a (2) + .... (15) 

7T \ 71 J 

The first order correction has been computed in || up to the finite part. Since it involves 
a second order pole in e and we are also interested in the finite part of its square, we need 
its expansion up to e 2 : 

« (1, ^(-0{^^c 2+£ (-I + Ic,)^(-H + fa)} + ^), 

(16) 
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where fi is the 'tHooft mass, and ( n = ((n) is Riemann's zeta function (£2 = tt 2 /6; 
£3 ~ 1.20206; £4 = 7r 4 /90). Furthermore, it is understood that (— l) e = exp(+ie7r). 

To this we add the second order correction: 
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where, as before, ni is the number of light quark flavors, ni = 5. 

Ultra-violet renormalization of the strong coupling constant is given by a® = Z a (Ma s )-a s , 
where Z a is related to the (5 function of Eq. (^) through 



— lnZ a (a s ) = — - 

oa s e-p{a s 



(18) 



5 Ratio of time-like to space-like form factor 



Both as a check and as an estimate on the magnitude of the corrections, we may consider 
the ratio of the time-like to the space-like form factor. It is free from infra-red singularities 
and contains the presumably most significant contributions stemming from the analytic 
continuation of the factor (^ 2 /(— q 2 )) e from space-like to time-like values of q 2 . 

As was shown in |2f| for the quark form factor in QCD |l6[ , this ratio can be derived up 
to 0(a 2 ) from the one- loop terms of the form factor by combining them with a known 
two-loop anomalous dimension. In the case of gg — > H this anomalous dimension is given 
by 9/4 times the one given in [25|. Following the derivation of [25] and setting /j, 2 
we arrive at the following expression: 
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(19) 



where we inserted n\ = 5 in the second line. The third line displays separately the LO, 
NLO and the NNLO contribution, as well as their sum, for a s (M^) = 0.112. Eq. (|l9|) 
fully agrees with a direct evaluation of the ratio using the expressions ([D]) , ( |l6| ) , and Ql7D 
for a. This provides a check on the terms cx a™/e 2n ~ fc (k = 0,1,2) of our result. The 
numerical value of the NLO correction in Eq. ( |l9| ) reflects the largeness of the full NLO 
terms as obtained in ||. The number for the NNLO corrections gives some hope towards 
a certain degree of convergence for the perturbative series of the full result for o~{gg — > H). 
Concluding this section, let us note that an interesting extension of this discussion could 
be the resummation of the leading terms along the lines of [p6[. 
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6 Conclusions and Outlook 



We used the recently introduced method of Q in order to calculate the NNLO virtual 
corrections to the production cross section of Higgs bosons in gluon fusion. The result is 
a gauge invariant component of the full cross section. The next step towards a complete 
answer for the NNLO rate is to integrate the squared amplitudes for the real radiation 
processes over the phase space. This is work in progress [^7). Finally, one has to convolute 
the full partonic cross section with the parton distribution functions. Their evaluation to 
the relevant order is therefore certainly a very important task. 
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